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ABSTRACT 

We present the results of the analysis of a large database of X-ray observations of 22 
galactic black-hole transients with the Rossi X-Ray timing explorer throughout its op- 
erative life for a total exposure time of ~12 Ms. We excluded persistent systems and the 
peculiar source GRS 1915+105, as well as the most recently discovered sources. The 
semi-automatic homogeneous analysis was aimed at the detection of high-frequency 
(100-1000 Hz) quasi-periodic oscillations (QPO), of which several cases were previ- 
ously reported in the literature. After taking into account the number of independent 
trials, we obtained 11 detections from two sources only: XTE J1550-564 and GRO 
J1655-40. For the former, the detected frequencies are clustered around 180 Hz and 
280 Hz, as previously found. For the latter, the previously-reported dichotomy 300-450 
Hz is found to be less sharp. We discuss our results in comparison with kHz QPO in 
neutron-star X-ray binaries and the prospects for future timing X-ray missions. 

Key words: accretion, accretion discs - black hole physics - relativistic processes - 
X-rays: binaries 



1 INTRODUCTION 

The Rossi X-Ray Timing Explorer (RXTE) mission has 
provided thousands of high-quality observations of black- 
hole transients (BHTs), which have profoundly changed our 
knowledge of the properties of accretion onto stellar-mass 
black holes (see e.g. Belloni 2010; Fender 2010; Belloni et al. 
2012). An important aspect has been the detection of Quasi- 
Periodic Oscillations (QPO) at frequencies higher than ~40 
Hz, which opened a new window onto fast phenomena in 
the frequency range expected from signals associated to Ke- 
plerian motion in the innermost regions of an accretion disk 
around a black hole. A small number of detections is avail- 
able in the literature from a a small number of sources: GRS 
1915+105 (Morgan et al. 1997; Strohmayer 2001a; Belloni 
et al. 2001; Remillard et al. 2002b; Belloni et al. 2006), GRO 
J1655-40 (Remillard et al. 1999; Strohmayer 2001b), XTE 
J1550-564 (Homan et al. 2001; Miller et al. 2001; Remillard 
et al. 2002a), H1743-322 (Homan et al. 2005; Remillard et 
al. 2006), XTE J1650-500 (Homan et al. 2003), 4U 1630-47 
(Klein- Wolt, Homan & van der Klis 2004), XTE J1859+226 
(Cui et al. 2000) and IGR J17091-3624 (Altamirano & Bel- 
loni 2012). The signals are very weak and some detections 
appear to have marginal statistical significance. Although 
these High-Frequency QPOs (hereafter HFQPOs) are very 
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important as they constitute a promising way to detecting 
effects of General Relativity in the strong field regime, no 
systematic work on the existing bulk of RXTE data is avail- 
able. In particular, it is important to compare quantitatively 
the existing results with those from kilohertz QPOs from 
neutron star X-ray binaries, where signals are stronger and 
the known phenomenology much better covered (see e.g. van 
der Klis 2006 for a review). The similarities between the tim- 
ing features of these two classes of sources are strong (see 
Wijnands et al. 1999; Psaltis et al. 1999; Belloni et al. 2002; 
van der Klis 2006) , but it is not clear how HFQPOs fit in this 
comparison. In this paper, we present a systematic analysis 
of a large sample of RXTE observations of BHTs, concen- 
trating on the detection of HFQPOs. All major sources are 
present, with the exception of GRS 1915+105, which will be 
the subject of a separate paper due to the large number of 
available observations and its peculiarities) and IGR J17091- 
3624, which was observed by RXTE after the analysis part of 
this work was completed. An a posteriori bulk analysis of the 
data suffers from problems related to number of trials, which 
means that the sensitivity to detections will be reduced. In 
addition, our procedure was aimed only at detecting one 
peak in each Power Density Spectrum (PDS), which means 
that double detections do not appear. Finally, different ob- 
servations were analyzed separately without any merging, 
although merging was the way some of these oscillations 
were found. However, such a (semi-automatic) procedure is 
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useful to assess what is available from RXTE data and to 
plan future observations with upcoming facilities like AS- 
TROSAT and LOFT. The comparison with previous works 
is in many cases difficult as the starting conditions for the 
analysis are markedly different. Motivated selection of ob- 
servations and/or merging, as well as the relaxation of the 
constraints imposed in our search (total number of available 
observations, frequency range considered, energy range for 
the accumulation, limits on the quality factor of the signal) 
can lead to more sensitive results. 



2 DATA AND ANALYSIS 

We selected all RXTE observations of known transient 
black-hole binaries available in the archive from the start 
of the mission until MJD 55601 (2011 February 9), con- 
centrating on the data from the Proportional Counter Ar- 
ray (PCA) instrument. We analyzed 22 sources, for a total 
of 7108 observations (see Tab. [l}. The bright source GRS 
1915+105 was not included for two reasons: its peculiarity 
and the fact that effectively during the RXTE lifetime it be- 
haved as a persistent source. The analysis of data from GRS 
1915+105 will be reported in a forthcoming paper. For the 
project presented here, it is important to establish a detailed 
analysis procedure in advance, in order to control and min- 
imize the probability of spurious detections. The following 
multi-step procedure was applied. 

(i) Two energy bands were considered for analysis: ab- 
solute channels 0-79 and 14-79 (corresponding to energies 
1.51-21.30 keV and 4.07-21.30 keV at the beginning of the 
mission and to energies 2.06-33.43 keV and 6.12-33.43 keV a 
the end of the mission, due to gain changes in the detectors). 
We will refer to these as the total and hard band respectively. 
For some observations, these exact bands were not available 
due to the data modes, in which case we chose the closest 
approximation to those boundaries. In a few cases, the data 
modes were such that we used the full channel range and/or 
the total and hard band resulted to be identical. The choice 
of the hard band was made in order not to deplete too much 
the number of counts and allow the examination of more 
observations (see next bullet). 

(ii) When analyzing a large number of observations, the 
issue of number of trials must be considered carefully. Of 
all observations in our sample, a sizable number were of a 
source at the very start or at the end of an outburst, where 
count rate is low and the sensitivity to detect quasi-periodic 
features is very low. We decided not to analyze observations 
where the data are not expected to yield a positive detec- 
tion for HFQPO. Rather than selecting a simple threshold 
in count rate to exclude the observations, we used the ex- 
pression to estimate the significance (in number of sigmas, 
Tier) of a broad feature in frequency domain 

1 2 / T eX p 

Ua ~ 2 rs (S + B) V Ai7 

(see van der Klis 1998) where r s is the fractional rms 
variability of the feature, S is the source net count rate, 
B is the corresponding background count rate, T exp is the 
exposure time and Au is the FWHM of the feature. We 



assumed a QPO with 6% fractional rms and 10 Hz FWHM 
and used the observed count rates and exposure time to 
estimate the expected sensitivity. We then excluded all hard 
or total observations which yielded n a < 3. 

With this procedure we retained 4205 total observations 
and 3106 hard observations, for a cumulative exposure time 
of 12.2 Ms and 10.0 Ms respectively (see Tab.[T|. We consider 
a separate observation (or dataset) the data corresponding 
to a RXTE observation ID. 

(iii) The large number of observations made it necessary 
to perform an automatic search in order to select candidates 
for HFQPOs. For each observation, the adopted procedure, 
applied both to the total and hard energy band, was the 
following: 

• We produced a PDS by dividing the dataset in in- 
tervals of 16s duration and averaging the corresponding 
PDS. The time resolution for all spectra was 4096 points 
per second, corresponding to a Nyquist frequency of ~2 
kHz. The PDS were normalized according to Leahy et al. 
(1983). The powers were rebinned logarithmically in such 
a way that each frequency bin was larger than the previ- 
ous one by ~2%. Uncertainties in the power values were 
estimated following van der Klis (1988). 

• Since the PDS of BHTs are rather complex and we 
were interested only in the high-frequency region, we con- 
sidered only the frequency range 100-1000 Hz. This means 
that our implicit definition of HFQPO is limited to cen- 
troid frequencies larger than 100 Hz. 

• We fitted the PDS with an automatic procedure with 
a model consisting of a power law (to account for the Pois- 
sonian noise component) and a Lorentzian. The slope of 
the power law was not fixed to (constant component) in 
order to fit any possible remaining tail of low-frequency 
features. We did the fit by first fixing the Lorentzian cen- 
troid frequency to all values between 100 and 1000 Hz in 
steps by 1 Hz. The FWHM of the Lorentzian was bound 
between 0.5 Hz and 1000 Hz. A fit with a free centroid 
was then made around the frequency corresponding to 
the minimum chi square. This resulted in a single best fit 
for each observation. Notice that this procedure is aimed 
at detecting a single HFQPO in each PDS and would not 
be able to detect multiple peaks. 

• From each fit, we estimated the significance of the 
detection by dividing the normalization (integral) of the 
Lorentzian by its la error on the negative side. Only 
features with this figure of merit larger than 3 were con- 
sidered. The surviving number of observations can be seen 
in the 3cr columns of Tab. [T] 

• Examining the best fit parameters, we were able to 
eliminate obvious meaningless fits. Those included fits 
where the FWHM had reached one of the two imposed 
limits (see above), making the feature impossibly narrow 
or broad, and those where the centroid frequency was close 
to the search limits. In addition, we discarded all features 
with a quality factor Q — vq/FWHM significantly less 
than 2, as we are interested only in peaked signals. We 
checked all these cases visually in order not to miss pecu- 
liar features. None were found. In this way, the numbers 
of HFQPO candidates dropped to those in the "good" 
columns in Tab. [l] 

• For all remaining PDS, we made a manual fit with 
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Table 1. Statistics of available data and obtained detections. Columns are: source name, total number of observations, number 
of observations analyzed in each band, total exposure considered in each band, number of observations with single try probability 
>3cr in each band, number of detections after manual analysis in each band and final number of detections. 



Source name 


N. obs. 


Valid T. 


Valid H. 


Exp. T (ks) 


Exp. H (ks) 


> 3<r T 


> 3<t H 


Good T 


Good H 


Final 


GX 339-4 


1353 


893 


517 


2190 


1657 


31 


34 


i 


11 


1 


4U 1630-47 


1002 


846 


701 


1934 


1672 


31 


27 


11 


7 


4 


GRO J1655-40 


596 


516 


491 


2295 


2214 


18 


24 


7 


8 


11 


H 1743-322 


504 


386 


300 


1153 


941 


19 


20 


6 


8 


1 


XTE J1550-564 


109 


319 


280 


860 


790 


16 


31 


15 


14 


16 


Swift J1753. 5-0127 


278 


240 


172 


671 


527 


14 


8 


3 


3 





XTE J1752-223 


208 


138 


79 


343 


250 


6 


6 





4 


2 


XTE J1650-500 


182 


101 


69 


219 


155 


17 


5 


1 


2 





Swift J1539.2-6227 


156 


61 


12 


153 


122 


1 


3 


1 


1 





XTE J1817-330 


155 


112 


79 


372 


256 


3 


2 


1 


1 





XTE J1859+226 


131 


115 


102 


314 


281 


6 


4 


2 


2 


1 


4U 1543-47 


104 


53 


38 


161 


125 


6 


3 





1 





XTE J1720-318 


100 


75 


19 


233 


101 


1 


1 





1 





4U 1957+115 


100 


97 


16 


529 


334 


1 














XTE J1118+480 


94 


62 


50 


171 


154 





3 





1 





MAXI J1659-152 


66 


63 


54 


148 


135 


1 


1 





1 





SLX J1746-331 


65 


17 


23 


125 


72 


2 


2 











XTE J1652-453 


57 


30 


1 


67 


2 


1 














XTE J1748-288 


24 


23 


21 


103 


99 





1 











SAX J1711. 6-3808 


17 


17 


12 


43 


30 

















GS 1354-644 


8 


7 


7 


52 


52 





1 











GRS 1737-31 


5 


1 


3 


40 


28 

















Total 


7108 


4205 


3106 


12177 


9996 


204 


176 


54 


65 


42 



the same model described above and we re-established the 
significance of the fit in the same way. In case we found a 
peak at a consistent frequency, both in the hard and total 
PDS, we retained the one with the highest significance. 
The final numbers of detections from this procedure are 
shown in the last column of Tab.Q] The total number was 
42. For all fits, the reduced chi square was close to unity. 

(iv) The resulting HFQPOs have a significance (esti- 
mated in the way described above) larger than 3cr for a single 
trial. We searched them in a number of PDS for each source 
and we considered a broad range of frequencies. Both these 
procedures translate in a number of trials. For each source, 
we considered how many PDS were searched (for the to- 
tal plus hard band) and how many independent frequencies 
were considered: for each detection we divided the 900 Hz 
interval used for the search by the FWHM of the detected 
QPO. Keeping these two effects in mind of course lowers the 
significance. 

We did this under the assumption that any frequency 
within the considered range would be accepted. If one wants 
to limit to certain frequencies, for instance considering only 
peaks around 180 Hz and 280 Hz for XTE J1550-564 as 
previously detected (Homan et al. 2001; Miller et al. 2001; 
Remillard et al. 2002a), the number of trials will change 
and so will the significance of the detections. For this rea- 
son, in Tab.[2]we show all 42 detections. In the following, we 
will restrict the discussion only to those with a high signifi- 
cance. We marked in boldface in the table the HFQPOs with 
a final chance probability less than 10 -2 . Notice that a num- 
ber of detections have a (rounded) final chance probability 
of 1, making them unlikely to be real. One has a rounded 



probability of 0, as it was so low that the spreadsheet we 
used could not approximate it. 

(v) Notice that our procedure can only detect one peak 
per PDS. In principle, two different features can be detected 
in the total and hard PDS, but this did not happen. No clear 
additional peaks were seen visually in the PDS containing a 
detection. 

In addition, we produced a Hardness-Intensity Diagram 
(HID) for each source, one point per observation, by accu- 
mulating net energy spectra from PCU2 (the best calibrated 
of the five units of the PCA) using standard procedures]^] 
The count rates were corrected for gain changes by divid- 
ing the spectra by a Crab spectrum (power law with photon 
index 2 and normalization 10 ph/cm~ 2 s _1 keV _1 at 1 keV) 
simulated with the detector response at the day of the ob- 
servation. Intensity is accumulated over PCA channels 8-49, 
corresponding to 3.2-19 keV at the end of the mission. Hard- 
ness is denned as H/S where 5* and H are the corrected rates 
in the channel bands 8-14 (3.2-6 keV) and 15-24 (6-10 keV) 
respectively. 



3 RESULTS 

In this section we examine the results for each source for 
which HFQPOs have been detected. We study in detail only 
detections at more than 3er significance after keeping number 
of trials into account. 



1 http://heasarc.nasa.gov/docs/xte/data_analysis.html 
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Table 2. List of detected HFQPOs. Columns are: observation date (MJD), exposure time, band in which 
the QPO was detected (Total/Hard), centroid frequency, FWHM and fractional rms of the QPO, single trial 
and final chance probabilities, source observed count rate (in PCU2) and X-ray hardness. 



MJD 


Exp. (s) 


B 


I/O (Hz) 


FWHM (Hz) %rms 


Po 


Pf 


Rate HR 


XTE J1550-564 



51076.000 
51101.607 
51108.076 
51115.281 
51241.802 
51242.507 
51245.354 
51247.979 
51255.158 
51258.497 
51259.253 
51291.184 
51664.409 
51664.637 
51665.406 
51668.829 



2960 
1584 
9872 
2048 
3104 
1504 
2752 
2816 
4080 
1104 
880 
1536 
2240 
2864 
2144 
2976 



T 
T 
T 
H 
H 
H 
T 
T 
H 
H 
H 
T 
H 
H 
H 
H 



180 

140. 
182 
272 
284 
280 
180. 
185. 
280 
276 
274 
134. 
274 
274 
264 
263 



70 
.27 
43 
07 
.00 
75 



- R +3.92 
•'° -4.37 

+ 1.71 
1.92 
+3.25 
3.20 
+3.42 
4.09 
-1.35 
-1.36 
+2.87 
2.49 
+2.02 
2.54 
9n +6.73 
zu -4.93 

84 -1.92 
n9 +6.64 
uz -6.20 
a _ +4.30 
8 ' -3.61 
»- +1.76 
'° -1.45 
+3.43 
3.47 
97 +2.90 



40 



-3.21 
h7.09 
-3.19 
1-9.41 
-9.96 



95.12 


+ 16.12 

-22.00 


1 


.15 


+0 
-0 


14 
.12 


16.51 


+ 19.85 
-4.84 


1 


.12 


+0 
-0 


59 
17 


77.88 


+ 11.65 
-9.77 


1 


.74 


+0 
-0 


20 
16 


30.59 


+ 16.30 
-8.95 


4 


.02 


+0 
-0 


78 
.47 


26.86 


+3.88 
-3.12 


3 


.66 


+0 
-0 


21 
17 


30.00 


+9.80 
-6.85 


3 


.78 


+0 
-0 


46 
32 


16.76 


+8.71 
-6.21 


1 


.08 


+0 
-0 


25 
.13 


46.04 


+ 15.69 
-12.13 


1 


.32 


+0 
-0 


30 
17 


24.37 


+8.53 
-6.28 


4 


.35 


+0 
-0 


50 
41 


36.77 


+23.00 
-12.36 


5 


.22 


+0 
-0 


01 
76 


26.00 


+ 14.57 
-8.68 


5 


.09 


+1 

-0 


00 
.67 


7.81 


+5.16 
-3.25 


6 


.84 


+1 
-1 


28 
09 


41.96 


+ 19.88 
-10.66 


6 


.66 


+1 

-0 


23 
64 


33.51 


+ 13.97 
-8.16 


5 


.42 


+0 
-0 


87 
51 


42.87 


+22.01 
-15.43 


6 


.50 


+1 
-1 


35 
02 




+30.55 


5 


.50 


+1 


56 




-17.58 




— 





2 


.87E-07 


1.62E-03 


15700 





68 


6 


64E-04 


1.00E+00 


3135 


0. 


.57 


1 


90E-08 


1.31E-04 


3636 





.53 


7. 


80E-06 


1.28E-01 


1785 


0. 


41 


8. 


77E-29 


0.00E+00 


4258 





.42 


2 


05E-09 


3.69E-05 


4107 





.43 


8. 


54E-06 


2.40E-01 


4718 


0. 


.54 


3 


30E-05 


3.21E-01 


4509 





.54 


7 


20E-08 


1.59E-03 


2329 





38 


2 


.91E-04 


9.86E-01 


1816 


0. 


35 


7 


53E-05 


7.90E-01 


1710 





39 


8. 


45E-04 


1.00E+00 


88 





.42 


1 


.11E-07 


1.42E-03 


1347 





29 


4. 


65E-08 


7.47E-04 


1419 


0. 


31 


6 


.87E-04 


1.00E+00 


1105 





.28 


1 


36E-04 


7.44E-01 


856 


0. 


26 



GRO J1655-40 



50296.311 


9056 


T 


273.72 ±1^1 


70.33 


50301.665 


6464 


T 


307.16 ±lrl 


68.46 


50311.391 
50317.441 
50324.380 


1856 
6176 
4896 


T 
H 
T 


285.56 tfhg 
443.16 tUl 
289.67 tt f 3 


63.78 
41.44 
98.78 


50330.254 
50335.913 


6144 
8160 


H 
H 


453.72 tfll 
446.35 till 


37.69 
38.74 


50383.565 


6192 


H 


442.34 till 


31.50 


50394.884 


5792 


T 


313.45 + m 


98.11 


53498.415 


13968 


T 


125.81 tlfo 


7.49 


53508.507 


7504 


T 


266.63 t 7 7 lo 


61.76 



+ 19.40 


1. 


.15 


+0 


15 


-10.92 


-0 


1)8 


+0.00 


0. 


.88 


+0 


21 


-38.07 




-0 


.11 


+42.48 


1. 


,14 


+0 


36 


-17.38 


-0 


16 


+20.27 
-13.66 


5. 


39 


+0 
-0 


81 
64 


+40.67 
-17.67 


1. 


.41 


+0 
-0 


30 
,11 


+ 16.46 


5. 


69 


+0 


89 


-12.16 


-0 


71 


+9.80 
-8.00 
+ 15.62 
-10.61 


6 


.03 


+0 
-0 


59 
54 


4 


.55 


+0 
-0 


70 
1,2 


+38.42 


1 


.17 


+0 


25 


-19.10 


-0 


.11 


+3.86 
-2.91 


0. 


65 


+0 
-0 


11 

.09 


+30.61 


0. 


.76 


+0 


19 


-13.60 


-0 


.09 



1 


15E-13 


1.483E-9 


7443 


0.67 


3 


.91E-05 


4.04E-01 


6414 


0.61 


1 


93E-04 


9.35E-01 


5450 


0.54 


1 


28E-05 


2.44E-01 


5210 


0.51 


1 


.59E-10 


1.46E-06 


7267 


0.65 


2 


.91E-05 


5.03E-01 


4964 


0.49 


9 


.01E-09 


2.11E-04 


4983 


0.50 


1 


.31E-04 


9.77E-01 


6092 


0.54 


3 


72E-08 


3.44E-04 


5497 


0.59 


2 


16E-04 


1.00E+00 


3190 


0.29 


2 


.67E-05 


3.24E-01 


10164 


0.57 



XTE J1859+226 



51485.875 



6064 



H 



269.44 tlfs 



38.55 tll fs 



4 91 + 1 ' 57 

-0.71 



2.70E-04 7.46E-01 



824 



0.21 



H 1743-322 



52803.517 
52805.425 
53207.602 
55422.024 



6416 
5232 
2608 
1440 



H 
H 
H 

T 



240.39 
232.38 



+3.99 
-4.26 
+ 1.41 
-1.38 

178.57 tl f 3 
208.36 tl'.ls 



33.11 
7.65 
6.01 
8.04 



+ 14.35 

-10.41 

+6.43 

-3.20 

+4.09 

-3.80 

+4.68 

-3.38 



3.01 
2.09 
6.84 
5.39 



+0.59 
-0.36 
+0.49 
-0.31 
+ 1.26 
-1.10 
+0.51 
-0.75 



1.66E-05 2.66E-01 1366 0.34 

3.13E-04 1.00E+00 1314 0.35 

9.04E-04 1.00E+00 316 0.14 

1.72E-04 1.00E+00 238 0.81 



GX 339-4 



52387.556 
52490.229 
52690.841 
54291.443 



11200 
816 
3408 
1264 



H 
H 

T 
T 



127.70 
278.67 
105.35 
208.44 



+4.16 
-5.57 
+3.47 
-3.09 
+ 1.77 
-1.50 
+2.81 
-2.68 



44.21 tllfs 

15-48 tl° 5 i 7 
+4.46 
-3.78 
+8.67 
-5.19 



8.87 
14.68 



+ 1.05 
-0.35 
1-1.45 
-1.01 
(-1.12 
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3.1 XTE J1550-564 

We selected sixteen features with a single-trial chance prob- 
ability lower than that corresponding to 3a: 1.35xl0 -3 (see 
Tab. [2}. When we took into account the number of trials 
and considered detections of significance higher than 99%, 
we were left with 7 cases. Notice that some detections cor- 
respond to a very low count rate, but we nevertheless in- 
cluded them in the sample for completeness. Most of them 
are at a frequency close to either 180 Hz or 280 Hz, with 
the exception of the ones with the highest chance proba- 
bility. In Figure [2] we show the HID of all observations of 
XTE J1550-564 with the detection of HFQPOs (depending 
on the chance probability). It is clear that all detections 
correspond to intermediate hardness and it is interesting to 
note that the observations lie on a diagonal line: the brighter 
the observed rate, the harder the energy spectrum, with the 
two detections at a lower frequency corresponding to softer 
spectra. The same applies to the gray points, with one ex- 
ception. An examination of the full PDS for these obser- 
vations indicates that the first one in Tab. [2] corresponds 
to the brightest observation, where there is a low-frequency 
QPO (LFQPO) which can be classified as type-A (although 
Remillard et al 2002a did not classified it at the time). The 
second observation shows a very clear type-B QPO, while all 
the others have a type-A QPO. Overall, all observations can 
be considered part of the soft-intermediate state (SIMS, see 
Belloni 2010). The centroid frequencies are clustered around 
two values: ~180 Hz and ~280 Hz, while their Q value is 
around 2 for the low- frequency peak and 6-10 for the high- 
frequency one. The fractional rms of the low-frequency peak 
is low, around 1% (total band), as opposed to that of the 
high-frequency peak which ranges from 3.6% to 6.7% (hard 
band). Notice that the two last observations belong to a dif- 
ferent outburst than the first five, but the QPO centroid 
frequency is remarkably similar to the others. 

There are 20 detections reported in Remillard et al. 
(2002a) from the 1998 outburst of XTE J1550-564. Ten of 
them appear also in our sample, whereas ten were not de- 
tected by our procedure. We checked them all visually and 
found that many of them resulted in very broad weak ex- 
cesses superimposed to the tail of a noise component, which 
led to non significant detections in our analysis. Remillard 
et al. (2002a) do not report any parameter other than the 
centroid frequency and Q values (from 1 to 15 from their Fig. 
7), and it is therefore difficult to make a deeper comparison. 
Remillard et al. (2002a) claim detections at 4<r level, but 
they give no indication as to how they estimated the signif- 
icance. 

Miller et al. (2001) report six detections of HFQPO 
from the 2000 outburst of XTE J1550-564. After taking into 
account a one-day shift in their Tab. 1, the first four of 
them correspond to our last four detections. The last two, 
at slightly lower frequency, are not in our final sample. An 
inspection of the PDS from those observations show that 
there is a feature, which our program detected at a signifi- 
cance slightly lower than 3c 
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Figure 1. PDS of the seven good detections of HFQPO for XTE 
J1550-564 with the corresponding best fit models. 
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Figure 2. Hardness-Intensity Diagram for XTE J1550-564. Each 
point corresponds to a PCA observation. Gray circles correspond 
to HFQPO detections with final chance probability larger than 
1%, black stars to those with smaller chance probability. 



3.2 GRO J1655-40 

This is the only other source of our sample with multiple 
overall detections: 4 out of 11 single-trial detections. The 
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Figure 4. Hardness-Intensity Diagram for GRO J1655-40 (see 
Fig.§. 



other observations from Strohmayer (2001b) appear in our 
sample, although for MJD 50311 the high-frequency peak 
found by Strohmayer appears slightly below our 3rr level. 



Figure 3. PDS of the four good detections of HFQPO for GRO 
J1655-40 with the corresponding best fit models. 



frequency of the QPO is also bimodal: between 273 Hz and 
313 Hz for the low-frequency peak (Q between 3 and 4, rms 
slightly above 1%) and 446 Hz for the high-frequency peak 
(Q=12 and 6% rms). From Fig. [4] it is evident that all de- 
tections correspond to the "anomalous state" at very high 
fluxes (see Belloni 2010). The PDS is also typical of that 
state. As in the case of XTE J1550-564, the HID points are 
distributed along a diagonal line. Overall, the four detec- 
tions with high significance after correction for number of 
trials show quite a range of frequencies: 274, 290, 313 and 
446 Hz. 

Detections of HFQPO in the 1996 outburst of GRO 
J1655-40 were reported by Remillard et al. (1999). They 
detect six features in single observations, in addition to the 
analysis of sets of combined observations, obviously not com- 
parable with our results. Four of their detections correspond 
to entries in our Tab. [2] one corresponds to our 443 Hz ob- 
servation at MJD 50317, but at a different frequency (294 
Hz); our detection is from the hard band, while theirs is 
from the total band. In the total band we do detect their 
signal, but just below 3a single trial. The same is for the ob- 
servation at MJD 50383. As in the case of XTE J1550-564, 
Remillard et al. (1999) give no details on the detections, nor 
on the procedure followed. Two of our detections do not ap- 
pear in their sample. The observation from MJD 50394 does 
not appear in their list, while that from MJD 50335 is a 
very significant peak at 446 Hz in our data, but in the hard 
band. Nothing can be seen in the total band, so our result is 
consistent with theirs. This last one is the observation where 
the ~450 Hz QPO was discovered (Strohmayer 2001b). All 



3.3 XTE J1859+226 

Only one feature appears in our sample, at a single trial 
probability of 2.7xl0 -4 , which translates in a very high fi- 
nal chance probability. This feature is very unlikely to be 
real. The source appears at intermediate hardness in the 
HID (see Fig. [5); the presence of a type-B QPO indicates 
that the source was in the SIMS. Cui et al. (2000) report 
three detections of HFQPOs. None of these QPO appears 
in our sample. Two of them appear to be below our > 3a 
single-trial threshold in their table (the significance is com- 
puted as the ratio between normalization and its lower-side 
error, while using fractional rms as done by Cui et al. (2000) 
increases the value by a factor of two because of the error 
propagation), while the third one is a very broad feature 
which is automatically discarded by our procedure. 



3.4 H 1743-322 

We detect four QPOs, with high final chance probability. 
From Fig. [5] we can see that in one case the source was 
in a very soft HSS, in one case the source was rather hard 
(the full PDS shows a type-C QPO, from which we conclude 
the source was in the HIMS), and in two cases the source 
showed intermediate hardness, and was in a similar location 
in the HID as the ones discussed for GRO J 1655-40 and 
XTE J1550-564. In these last two cases the full PDS shows 
a strong type-B QPO, and therefore the source was in the 
SIMS. HFQPOs were reported from H 1743-322 by Homan 
et al. (2005) at 163 and 240 Hz from the observations of MJD 
52787. They do not appear in our sample, but accumulation 
of more than one observation is needed to obtain a detection, 
so our results are not inconsistent with theirs. The same 
applies to the detection of Remillard et al. (2006) . 
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Figure 5. Hardness-Intensity Diagram for the remaining sources for which no HFQPO with chance probability < 1% was found after 
accounting for the number of trials. Observations with HFQPOs with chance probability > 1% are marked with thick circles and crosses. 



3.5 GX 339-4 and XTE J1752-223 

For GX 339-4, the fur detections form a square in the HID 
(see Fig.[5|. The two hard ones are obviously LHS. The two 
soft ones are HSS. For XTE J1752-223 the two detections 
are at the same high value of hardness, in the LHS. In 
these two case the chance probability was, respectively, 0.23 
and 0.035. These are most likely statistical fluctuations and 
therefore we do not discuss them further. 



3.6 4U 1630-47 

None of the four peaks in Tab. [2] appear to be statistically 
significant after the (large) number of trials is taken into 
account. Two peaks are detected when the source flux was 
low: in one case the source was in the LHS, while in the 
other case the full PDS shows very little noise, making it 
difficult to classify the state. The remaining two observa- 
tions lie along a diagonal branch in the HID and also show 
low noise level. High-frequency features in this source were 
reported by Klein- Wolt et al. (2004). Those are very broad 
features which depend crucially on the subtraction of Pois- 
sonian noise. However, all of them were obtained after aver- 
aging observations corresponding to different intervals in the 
1998 outburst and therefore cannot appear in our sample. 



3.7 XTE J1650-500 

We did not obtain detections for this source. Homan et al. 
(2003) report high-frequency features, some of which very 
broad like in the case of 4U 1630-47. Once again, consid- 
erable averaging of more observations has been performed, 
which explains the absence of detections in our sample. 



4 DISCUSSION 

We analyzed 7108 RXTE/PC A observations of 22 black-hole 
transients searching for HFQPO. From the whole sample, af- 
ter taking into account the number of trials due to the many 
observations for each source and the number of independent 
frequencies sampled, we obtain only 11 detections with a 
chance probability less than 1%. These detections belong 
to two sources : XTE J1550-564 (7) and GRO J1655-40 (4). 
The small number of detections (including single-trial detec- 
tions at significance higher than 3<r there are only 42), con- 
firms that these are very elusive signals, as already indicated 
by the small number of detections reported in the literature. 
A few additional detections were obtained by other authors 
by averaging observations according to hardness, but those 
cases will not be discussed here. From our detections, a num- 
ber of aspects can be discussed. 

All detections correspond to PDS from very specific 
states. Six out of seven detections for XTE J1550-564 are 
from the SIMS, while the last one corresponds to the high- 
soft state (HSS). One observation has a clear type-B QPO, 
three have a type-A QPO and one corresponds to the bright 
flare in 1998, where a LFQPO can be tentatively identified 
with type-A (see Remillard et al. 2002a). No type-C QPO 
is detected simultaneous to one of the HFQPO. The first 
five HFQPO are from the 1998 outburst, the last two from 
the 2000 outburst (see Rodriguez et al. 2004). The four de- 
tections from GRO J 1655-40 correspond to the "anomalous 
state" described in Belloni (2010). Again, no type-C QPO 
can be identified in those four PDS. These detections come 
all from the 1996 outburst (Remillard et al. 1999). From the 
HIDs in Figs. [2] and [3] we see that the points corresponding 
to the HFQPO detections in both sources are distributed 
diagonally. However, it is difficult to relate this to precise 
spectral characteristics. 
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For XTE J1550-564, the centroid frequencies of the de- 
tected QPOs appear to cluster around two values: ~ 180 
Hz and ~280 Hz. Interestingly, the low-frequency detections 
come from the total energy band, while the high-frequency 
detections come from the hard band, indicating that the 
high-frequency oscillation has a harder spectrum. Moreover, 
inspection of Fig. [3] shows that the low-frequency peaks are 
clearly broader than the high-frequency ones. 

For GRO J1655-401, the dichotomy 300-450 Hz seems 
to be less sharp. Strohmayer (2001) detected the ~450 Hz 
feature simultaneous with the 300 Hz one when averaging 
observations, showing that the two peaks correspond to dif- 
ferent time scales. We see also frequencies (~270 and ~290 
Hz) significantly lower than 300 Hz. Also in this source the 
high-frequency detection is significantly narrower than those 
at lower frequencies. 

The fact that the few available detections appear at 
more or less the same frequencies for XTE J1550-565, even 
for different outbursts, and that the upper peak is nar- 
rower than the lower make a comparison with kHz QPO in 
neutron-star binaries difficult. There, the centroid frequen- 
cies are observed to vary in a way consistent with a random 
walk (see Belloni et al. 2005, 2007) and the upper peak is 
broader than the lower one (see e.g. Barret et al. 2005). 

Two sources which are known to show HFQPOs were 
excluded from the analysis. GRS 1915+105 is a very pe- 
culiar system with many observations and many HFQPO 
detections. We will present our results on this source in a 
forthcoming paper. IGR J17091-3624 is a new "twin" of GRS 
1915+105, in the sense that it is the only other known source 
to display its unusual variability. Results were already pub- 
lished in Altamirano et al. (2012). 

It is interesting to compare our results with those of kHz 
QPOs in neutron-star LMXBs (see van der Klis 2006 for a 
review). The neutron-star system 4U 1636-53 was observed 
regularly during the last few years of the RXTE lifetime, 
yielding a solid database from which to extract statistical 
information. In this source, ~40% of the observations led 
to the detection of at least one high-frequency feature (see 
Sanna et al. 2012) as opposed to 0.15% here. 

Revnivtsev & Sunyaev (2000) showed that there is a 
systematic difference in power level between black-hole and 
neutron-star sources in the hard state, attributing it to 
the presence of a radiation-dominated spreading layer in 
the latter, the so-called boundary layer. As already sug- 
gested by these authors, the presence/absence of strong 
high-frequency QPOs can be part of the same picture. Fol- 
lowing this line, Gilfanov et al. (2003) and Gilfanov & 
Revnivtsev (2005) analysed low-frequency and kHz QPOs in 
neutron-star binaries and concluded that the oscillations are 
caused by variations in the luminosity of the boundary layer. 
In this framework, it is the diminishing relative importance 
of the boundary-layer emission which weakens kHz QPO as 
they move to very high frequencies (see Mendez 2006 for 
additional discussion). Moreover, our HFQPOs are detected 
in intermediate states. Belloni et al. (2007) and Sanna et 
al. (2012) showed that in the neutron-star LMXB 4U 1636- 
53 kHz QPOs are not found in the observations when the 
source is the hardest, at intermediate hardness mostly the 
upper kHz QPO is detected and in the softest observations 
mostly the lower kHz QPO is present. It is difficult to make a 
comparison with our case: that no kHz QPOs are detected in 
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Figure 6. Frequency distributions for all detections of upper (top 
panel) and lower (bottom panel) kHz QPOs in the NS LMXB 4U 
1636-53 (Sanna et al. 2012) (gray histograms). The black his- 
tograms include only detections with rms >15% for the upper 
and >10% for the lower. 



the hard state is consistent with the black-hole case, while a 
detailed comparison between soft states remains to be done. 

The detection of only a few specific frequencies in the 
case of black-hole HFQPOs is something that does not ap- 
ply to kHz QPOs in neutron-star binaries. The frequencies 
of kHz QPO are distributed along a very broad range (see 
van der Klis 2006, Mendez 2006, Belloni et al. 2007, Barret 
et al. 2005, 2007). KHz QPOs in neutron stars are usu- 
ally strong, more easily detected than HFQPO, which then 
leads to a much larger number of detections. Moreover, at 
most frequencies the lower kHz peak is narrower and more 
significant than the upper one. The question is whether the 
frequencies of the oscillations can be understood in terms of 
the same physical mechanisms in the two classes of sources. 
In order to check what frequencies would be observed if this 
was the case, in Fig. [5]we plot the distribution of upper and 
lower kHz QPOs in the full RXTE database of 4U 1636-53 
(data from Sanna et al. 2012). In both panels, the black his- 
togram contains only the strongest detected peaks (we chose 
an rms threshold of 15% and 10% for the upper and lower 
peaks respectively in order to select only a handful of detec- 
tions). A comparison with the HFQPOs of XTE J1550-564 
is interesting. There, the lower-frequency HFQPO is weaker 
and broader than the higher-frequency one. Unfortunately, 
the only reported case of a double-detection in XTE J1550- 
564, obtained through averaging of many observations, does 
not allow a good determination of the quality factor of the 
two peaks (Miller et al. 2005). The black histograms show 
that on average the strongest upper kHz QPOs in 4U 1636- 
53 appear at a lower frequency than the strongest lower kHz 
QPOs. This is consistent with the spread of the frequencies 
we find in XTE J1550-564. Clearly, the limited number of 
detections of XTE J1550-564 does not allow a more precise 
comparison. 

The width of the QPO peaks is another interesting 
quantity of a comparison with neutron-star binaries. In our 



HFQPO in black-hole binaries 9 



detections, the Q factor of the upper-frequency peak is 
higher than that of the lower- frequency one. In neutron- 
star binaries it is generally the opposite (see e.g. Barret et 
al. 2007). This is a clear difference between the neutron-star 
and black-hole cases. If the process originating the oscilla- 
tions is the same of the two classes of courses, the presence 
of the boundary layer is likely to be the reason for this dif- 
ference, yielding another observational constraint for theo- 
retical modeling (see Mendez 2006). 

Overall, we detected eleven significant features after an- 
alyzing around 7000 observations, resulting in a detection 
efficiency of ~ 10~ 3 . Of course, since the HFQPOs are de- 
tected only in intermediate states and these states are the 
short-lived ones, observing with new instruments and con- 
centrating on the most likely periods during outbursts can 
increase the probability of detection. In other words, to de- 
tect HFQPOs one does not need a good coverage of a full 
outburst, but only of the intermediate states. The LAXPC 
instrument on board the upcoming ASTROSAT mission 
(Agrawal 1996) will increase the effective area of the PCA 
above 10 keV, where the HFQPOs are stronger (see Morgan 
et al. 2001), but in order to be effective ASTROSAT will 
have to concentrate observations during bright intermediate 
states. The LOFT mission (Feroci et al. 2010), under evalu- 
ation by ESA at the time of writing, is planned with a very 
large effective area (around 12 m 2 ) and will allow sampling 
of these weak signals in all states. 
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